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Majorana fermions are long- sought exotic particles that are their own antiparticles. Here we propose to 
utilize superconducting circuits to construct two superconducting-qubit arrays where Majorana modes can 
occur. A so-called Majorana qubit is encoded by using the unpaired Majorana modes, which emerge at the 
left and right ends of the chain in the Majorana- fermion representation. We also show this Majorana qubit 
in the spin representation and its advantage, over a single superconducting qubit, regarding quantum 
coherence. Moreover, we propose to use four superconducting qubits as the smallest system to demonstrate 
the braiding of Majorana modes and show how the states before and after braiding Majoranas can be 
discriminated. 

Majorana fermions are particles that are their own antiparticles. These long-sought particles have recently 
received considerable interest (see, e.g., Refs. 1-11). It has been recognized 12-14 that a relatively easy-to- 
engineer system — one-dimensional (ID) semiconducting wires on an s-wave superconductor — can 
realize a nontrivial topological state supporting Majorana fermions. This state is characteristic of ID topological 
superconductors 5 , in which Majorana modes can occur without requiring the presence of vortices in the system. 
The recent experimental observation 15 of a zero-bias peak in the differential conductance of a semiconductor 
nanowire coupled to a superconductor suggested the possible existence of Majorana fermions. Moreover, it was 
proposed 16 to use tunable ID semiconducting wire networks on an s-wave superconductor to demonstrate the 
non-Abelian statistics of Majorana fermions, because the Majoranas in the semiconducting wires can also behave 
like vortices in a p + ip superconductor 3,4 . In addition, it was also recognized 1718 that when a Jordan-Wigner 
transformation is performed, a ID quantum Ising model is equivalent to a ID topological superconductor, and 
Majorana modes can also occur therein. Nevertheless, less attention has been paid to this quantum Ising model 
than to ID topological superconductors because it was often regarded as a toy model. 

In this paper, we propose to realize such a toy model by using experimentally accessible superconducting-qubit 
arrays. Importantly, superconducting qubits can behave as controllable artificial atoms and tunable interqubit 
couplings are also achievable (see, e.g., Ref. 19-21). For instance, the tunable coupling between flux qubits was 
experimentally demonstrated in Refs. 22-24. For a finite superconducting-qubit array, when the interqubit 
couplings are tuned to be nonzero and other parameters of the qubits are tuned to be zero, there are two unpaired 
Majorana modes, which emerge at the left and right ends of the chain in the Majorana-fermion representation. 
We use these two Majorana modes to encode a qubit which is here called the Majorana qubit. Also, we express this 
Majorana qubit in the spin representation and show its advantage, over a single superconducting qubit, regarding 
quantum coherence. Moreover, the advantages of superconducting qubits in controllability make it possible to 
construct a tunable ID quantum Ising model on wire networks, similar to the semiconducting wire networks in 
Ref. 16, to demonstrate the non-Abelian statistics of Majorana modes. We propose to use four superconducting 
qubits as the smallest circuit to demonstrate the braiding of Majorana modes, and show how the states before and 
after braiding Majoranas can be discriminated. This should provide an experimentally realizable, relatively simple 
setup to manipulate and probe Majorana fermions. Thus, our proposal could allow the quantum simulation 25 or 
emulation of Majorana fermions. 

Results 

Majorana modes in superconducting circuits. We construct two types of superconducting-qubit arrays (see 
Figure 1), which can exhibit Majorana modes. 
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Figure 1 | Two arrays of superconducting qubits. (a) Charge-qubit array: 
Nearest-neighbor charge qubits Q n and Q n+ \ are coupled by a large 
Josephson junction with coupling energy E Jc (shown as a crossed 
rectangle), (b) Flux-qubit array: Nearest-neighbor flux qubits are coupled 
by a coupler consisting of a flux-biased loop that is interrupted by two large 
Josephson junctions (each with coupling energy E Jc ) and a small Josephson 
junction with coupling energy PE Jcy where 0 < /?<^Cl. In (a) and (b), is 
the flux applied to each qubit loop, (c) Main components of a charge qubit, 
where a superconducting island (denoted as a solid circle) is connected to 
two Josephson junctions (each with coupling energy Ej 0 <^Ej c and 
capacitance Cj) and biased by a voltage V g through a gate capacitance 
Cg<^Cj. (d) Main components of a flux qubit, where two Josephson 
junctions, each with coupling energy Ej<^E Jc , connects a symmetric dc 
SQUID biased by a flux 0 5 . 

(1) Charge-qubit array. For the array of charge qubits shown in 
Figure 1(a), every pair of nearest-neighbor qubits are coupled by a 
large Josephson junction acting as an effective inductance. The non- 
nearest-neighbor qubits can also be coupled via these large Josephson 
junctions, but the interactions are negligibly small. Here we assume 
that all charge qubits are identical and that all large junctions are 
equal to each other. When leading terms are considered, the 
Hamiltonian of this charge-qubit array can be written as 

tf=X>« + i-X><+K), a) 



with fi = ~E c h (l — CgVg/ e )> v = Ejo cos(nQ> q /Q>o), and the interqubit 



coupling is given by 2< 



t = Lr 



TlEjQ 

#0 



.2 f n % 

O 0 



(2) 



Here £ c h (~e 2 /C/) ^>Ej in the charging regime considered here and 
Lj = <S> 0 /2tiI c , with I c = 2nE Jc /(D 0 and <X> 0 being the flux quantum. The 
eigenstates of the Pauli operator a z n are the charge states |0„) and 1 1„), 
corresponding to zero and one extra Cooper pair in the supercon- 
ducting island of the nth qubit. The Hamiltonian (1) provides an 
analog to the ID quantum Ising model. 



We now consider the case with the fluxes in all charge-qubit loops 

being tuned to <S> q = - so that v = 0, and the interqubit couplings 

reach the maximum t = L^tlEjqI^q) 2 . Using the Jordan-Wigner 
transformation 1718 : 



n-l ., , 

a n = G~ IT a\ = a\ W o 

m—l ... - i 



n-l 

n 

m= 1 



(3) 



,±_ 1 



where = - (a x n ± ia y n ), one can cast equation (1), in the case of 
%=^$>o,to 

N-l N 

H= ^ t (a n -al)^a n+1 +al +1 j-^^(2ala n -l), (4) 

n= 1 n— 1 

where the Dirac fermions obey the anticommutation relation 
|a n ,a^| =S nn '. Introducing Majorana fermions: 

(5) 

one can rewrite the Hamiltonian (4) as 



y^ = a\ + a n , y B n = i(a\-a n ), 



H=iJ2t Y y n+1 -ij2M 



(6) 



where yf =f n and {y* 9 y%} =2dxx>d nn >. Obviously, (}£) 2 = 1, 
which is different from the Dirac fermion. 

(2) Flux-qubit array. Figure 1(b) shows an array of flux qubits. Here 
the small junction in the ordinary flux qubit is replaced by a sym- 
metric dc SQUID to increase the tunability of the qubit. Also, a 
coupler consisting of three Josephson junctions is used to produce 
a controllable interqubit coupling between nearest-neighbor flux 
qubits. We assume that the parameters are the same for all qubits 
and also for all couplers. Moreover, the plasma frequency of the 
coupler is much higher than the related qubit energy, so as to keep 
the coupler in the ground state 27 . When the leading terms are 
included, the Hamiltonian of the flux-qubit array can be written as 



H=^K< +1 -E( vff «+^)- 



(7) 



Here v = I p <$} 0 ^- —f J , with I p being the persistent current of the flux 

qubit and / = <3> q /<3> 0 + f s /2, where/, = <3> s /<3> 0 , with ® s being the 
magnetic flux applied in the SQUID loop [see Figure 1(d)]. The 
eigenstates of the Pauli operator a z n are the clockwise and anti- 
clockwise persistent-current states of the nth qubit. The symmetric 
SQUID provides an effective Josephson junction with coupling 
energy aEj, where a = cos(nf s ). The exact expression of fi in equation 
(7) cannot be obtained, but it depends on a; numerical results 28 and 
approximate analytical calculations 29 showed that /n = 0 when a = 1. 
The interqubit coupling strength reads 27 



t = 



PE Ic cos(2nf c -(l) c ) 
l+2jScos(27i/ c -^ c ) ' 



(8) 



where/ c = O c /O 0 is the reduced flux applied to the coupler, and <j> c = 
2/3 sin(27r/ c )/[l + 2/5 cos(2rc/ c )], with ft being the ratio of the 
Josephson couplings between the smaller and larger junctions in 
the coupler [see Figure 1(b)]. 

We study the case with/ = ^ for all flux qubits, so as to have v = 0. 

The Hamiltonian of the system also becomes equation (4) when 
applying the Jordan-Wigner transformation: 
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n — 1 n — 1 




where = - (<J z n ±i<J y n ). Finally, the Hamiltonian is described by 

equation (6) when introducing Majorana fermions in equation (5). 
Therefore, the resulting Hamiltonians in terms of Majorana fermions 
are the same for both charge- and flux-qubit arrays. 

For N — > °°, we can obtain the energy bands of the periodic chain 
by performing a Fourier transform on Hamiltonian (6): 

v v k 

where y x _ k = y^ and X = A,B. The resulting Hamiltonian in recip- 
rocal space reads 

»-?(* *)U oo 

with D(k) = te ik + [i. The energy spectrum shows the particle-hole 
symmetric dispersion 

E(k)=±\D(k)\=±\te? k + /i\, (12) 

which consists of two bands. As examples, we present in Figure 2 the 
particle-hole symmetric dispersion for r = /n/t= 0.5 and 1, respect- 
ively. It is clear that when r = 1, the gap of the two bands closes at 
certain values of the wave vector k. 

For a finite chain, when \r\ < 1, there are two degenerate edge 
modes with zero energy (i.e., in the middle of the energy gap). These 
two edge modes can be represented by 

Q = ciyf + c 2 yf+ • • • + c 2N -i)j I + c 2N )%, (13) 
with the coefficients determined by 

fic 2n -i + tc 2n + i = 0, tc 2n -2 + fic 2n = 0, (14) 

where n = 1,2, . . ., N, and the initial condition is c 0 = 0 for the left- 
end edge state and c 2 n + i — 0 for the right-end edge state. It can be 
derived that the left- and right-end edge modes are given, respect- 
ively, by 

QL = C[yA-ryj + r 2 yA- ■■■ +(-r) N -V N ], 



Q R = C[(-r) N -y i + ---+ r y N _ 2 -ry B N _ 1 +y B N }, (15) 

where the normalization factor is C = (^~^_* r 2 ^ . In Ref. 30, 

the left- and right-end edge modes were also studied in a flux-qubit 
array, but the interqubit coupling was not tunable. 

In particular, when fi = 0, the Hamiltonian is reduced to 

H=-Y / t(2dld n -l), (16) 

n=l 

where d n = ^ (yjj +1 + is a Dirac fermion composed of two 
Majoranas at adjoining sites. The edge modes become two unpaired 
Majorana fermions: Qi = yf , and Qr = y^, which emerge at the left 
and right ends of the chain as local modes in the Majorana-fermion 
representation. However, as shown in the following subsection, these 
two Majorana modes become non-local in the spin representation. 
Moreover, these two degenerate modes do not appear in the 
Hamiltonian because they have zero energy 5,31 . Now define \F) to 
be the state in which all eigenstates of the system with E < 0 are 
occupied and those with E > 0 are empty. When the edge modes are 
occupied, |^) =yf \F) and ^r) =y^\F) are two degenerate ground 
states of the system. These two Majorana modes can be used to 
represent the basis states of a qubit called here the Majorana qubit: 

|0}=d end |F>, |1> = 4J0>, (17) 

where 

*nd=^(yf + r$) (18) 

is a non-local Dirac fermion, and d end \0) = 0. A similar Majorana 
qubit was also proposed in quantum wires with spin-orbit interac- 
tions (see, e.g., Ref. 16). Such a qubit had initially been thought of as 
being fully topologically protected, but recent studies showed that it 
could also suffer from decoherence caused by either coupling to the 
solid-state environment (see, e.g., Ref. 32) or strong renormalization 
by interactions (see, e.g., Ref. 33-35), which was often neglected. 

The Majorana qubit in the spin representation. Below we derive 
the two basis states |0) = d end \F) and |l)=dj nd |0) in the spin 
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representation and discuss issues regarding the quantum coherence 
of this Majorana qubit. 

(1) Charge- qubit array. When t > 0 and fi = v = 0 in equation (1), 
the state \F) can be written, in the spin representation, as 



\F) = 



V2\- 
V2\- 



if N = even, 
if iV = odd, 



where N is the number of charge qubits in the array and 
1 fl\ , 4 1 / 1 



y/2\l 



y/2\-l 



(19) 



(20) 



are the two eigenstates of o x with eigenvalues 1 and — 1, respectively. 

In the spin representation, y\ = o\ and y^ = iP c (^, where the 
string operator 



.= n * 



(21) 



is the parity operator associated with the Z 2 symmetry of the system. 
Obviously, y^ becomes a non-local chain operator in the spin rep- 
resentation. The two degenerate ground states |*F L ) and |^) can be 
written as 



and 



\^L)=7t\F)=at\F) 
\V R )=f N \F)=iP c o* N \F) 



V2\< 



-) if N = even, 
>) ifiV = odd, 



>) if iV = even, 
if iV = odd. 



(22) 



(23) 



It is clear that CF^a) = 0. The two basis states |0) = d end \F) and 
|1) = dl nd \0) of the Majorana qubit are given, respectively, by 



\o) = ^ + iy B N )\F) 
i ,. 



>)-\< 



>)) if N = even, (24) 
-)) if2V = odd, 



and 

|i)= 



;(yf-*^)|o> = zK+Pc^)|o) 



1 
1 

V2 



(I- 



>+ 



>)) if N = even, (25) 
-)) ifN = odd. 



Note that |l) = dj nd |0) can also be written as |1) = d\ nd d Qn( \\F) = 
^ (l + /yfy^)|F), which is identical to equation (25). These two 

basis states of the Majorana qubit are also two degenerate ground 
states of the system. Moreover, these ground states have well- 
defined parities because P c |0) = |0) and P c |l) = — 11) if N = even, 
and because P c \0) = -|0> and P c \l) = |l) if N = odd. In Ref. 36, 
similar Majorana modes in spin-chain networks were also used to 
encode a qubit. 



When the externally- tunable parameters such as gate voltages and 
applied fluxes are identified at each qubit, equation (1) can be rewrit- 
ten as 

N-l N 

h= wi «+i - (»« < +v - <)' ( 26 ) 



where 



cX 



.(«) 



, v n = £ /0 cos 



o 0 



and the interqubit coupling is given by 

_ / 7L£}/ 0 " 2 
ln,n+\ — L J 



(„.„+i=L, I—— sin sin — 4- 

«D 0 J \ 0>o / \ <D 0 



, (27) 



(28) 



EN 
Vn&t ^ 0, this longitudinal term will lift 
n—\ n 

the state degeneracy of the system. However, in our designed circuits, 
this can be avoided because we can have v n = 0 for each qubit by 

tuning the external flux to <D^ = - <D 0 . Also, we can have \i n = 0 by 

tuning the gate voltage to vj^ =e/C g , so as to achieve unpaired 
Majorana modes emerging at the two ends of the charge-qubit array 
in the Majorana-fermion representation. 

With regard to the quantum coherence of the Majorana qubit, 
there are three types of local perturbations that we should consider: 
(i) 3}i n o z n , (ii) 3v n o*, and (iii) 3t n ^ n+ \ o x n a x n+l . The charge perturba- 
tion 3/n n o z n can be explicitly written as 



where 



3fin= ~ \2e) d( ^ n ' 



3Q n = C g 3V^ + 3Q^\ 



(29) 



(30) 



with the term C g dvj^ arising from the gate- voltage fluctuations and 

(n) 

3Q b being due to the back-ground charge fluctuations (e.g., the two- 
level fluctuators). As shown in equations (27) and (28), the para- 
meters v n and t nt „+i contain both the Josephson coupling E J0 and 
the flux <D^\ Therefore, the local perturbations 3v n o x n and 
St njn+ i o x n G x n+l can be contributed by both the critical -current 37 
and flux fluctuations. 

The local perturbation 3/n n a z n can only tend to drive the ground 
state (i.e., the Majorana-qubit state) |0) (| 1)) to an excited state, which 
has an energy level higher than the ground state. This is owing to the 
protection of the Majorana-mode states ^j) and against the 
local perturbation 3fi n o z n , because this perturbation cannot produce 
a state transition (relaxation) between and |^a). Actually, the 
local perturbation 3}i n o z n tends to drive ^j) (1^)) to an excited 
state with an energy difference A from the ground state, where A = 4t 
for 1 < n < N, and A = 2t for n = 1 and N. Nevertheless, such a state 
transition is not permitted for a small perturbation 3\i n o z n . Thus, the 
local perturbation 3fi n o z n (i.e., the charge fluctuations) will not pro- 
duce decoherence to the Majorana qubit. This is a distinct advantage 
of the Majorana qubit over a single charge qubit in which the charge 
fluctuations dominate. 

The environmentally-induced decoherence in a near-critical ID 
system of AT^>1 coupled qubits was studied in Ref. 38, where a model 
Hamiltonian analogous to equation (1) with v = 0 was used and only 
the local magnetic-field fluctuations (i.e., the local perturbation 
3fi n a z n in our model) were considered. It was found that the require- 
ment of preserving the qubits' entanglement over a certain idling 
time between consecutive gates can be better fulfilled away from 
criticality, i.e., when r(=fi/t) tM. In our study, we consider the case 
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with r = 0, which is away from the criticality, and the two Majorana- 
qubit states |0) and |1) are entangled states of multiple qubits [see 
equations (24) and (25)]. Indeed, as discussed above, these multi- 
qubit entangled states are robust against the local perturbation 

3fi n a z n . 

As for the local perturbations 3v n o x n and 3t n ^ n+ \ o x n a x n+v they 
should randomly shift the energy levels of the states and 
I^r), causing pure dephasing to these Majorana-mode states. 
However, while dt^n+i o x n a x n+l yields pure dephasing to the 
Majorana-qubit states |0) and 1 1), the local perturbation 3v n a x pro- 
duces relaxation to these Majorana-qubit states. In a circuit com- 
posed of inductively- coupled charge qubits, the interqubit coupling 
is usually much smaller than the Josephson coupling energy E J0 , so 
the coupler perturbation 3t n ^ n+ \ o x n o x n+l should be weaker than 
3v n o x n . As shown above, the Majorana qubit is robust again the 
charge noise. Now, the dominant noise in the Majorana qubit is 
due to the perturbation 3v n o x n involving both critical- current 
and flux fluctuations. In order to have a longer decoherence time, 
a single charge qubit usually works at the optimal (i.e., degeneracy) 
point N g = eV g /2e =1/2. When this single charge qubit is slightly 
away from the optimal charge degeneracy point, the decoherence 
time becomes drastically short because of its strong sensitivity to 
the charge noise. Nevertheless, the Majorana qubit consisting of a 
charge-qubit array is robust against the charge noise. Then, its 
quantum coherence is still preserved even if each charge qubit is 
randomly shifted away from the optimal charge degeneracy point. 
This is also one of the advantages of the Majorana qubit over a single 
charge qubit. 

In Ref. 39, an inhomogeneous spin ladder was proposed to study 
the robustness of the Majorana modes. This spin model is an 
inhomogeneous ladder version of the Kitaev honeycomb model 40 . 
Similar to the ID quantum Ising model, the zero- energy Majorana 
modes of the inhomogeneous spin ladder are also localized in 
the fermionic representation and emerge at either the two ends of 
the ladder or the boundary between sections in different topological 
phases 39 . As shown above, in the quantum Ising model described 
by equation (1), the topological ground- state degeneracy is 
robust against the local perturbation 3fi n o z n , but can be lifted by 
the local perturbation 3v n a x . In the inhomogeneous spin ladder, 
the topological ground-state degeneracy cannot be fully lifted 
by inhomogeneous magnetic fields purely along the x, y or z dir- 
ection 39 . This is the advantage of the inhomogeneous spin 
ladder. However, as further shown in Ref. 39, the topological 
ground-state degeneracy of the inhomogeneous spin ladder can be 
lifted by local two-body terms. In the ID quantum Ising model 
in equation (1), the two-body (i.e., coupler) perturbation 
3t n ^ n+ i o x n a x n+l can also lift the topological ground-state degeneracy, 
but compared with the local perturbation 3/n n a z n , the local perturba- 
tion 3v n a x n is much weaker and the coupler perturbation 
St n>n +i o x n o x n+l is even weaker in the ID quantum Ising model rea- 
lized using a charge-qubit array. 

As a variation of the charge qubit, the transmon qubit was also 
often used in superconducting quantum circuits 41 . In this qubit, the 
perturbation 3v n g x , which can be due to the fluctuations of flux, 
cavity photons and critical current, is more important than the per- 
turbation 3fi n o z n arising from the charge noise. In the Majorana 
qubit with charge qubits replaced by transmons, the perturbation 
3v n <j x becomes more important, but the advantage of the 
Majorana qubit regarding the insensitivity to the charge noise will 
still remain. Therefore, the quantum coherence of the Majorana 
qubit is preserved even if each transmon shifts randomly away from 
the optimal charge degeneracy point. 

(2) Flux-qubit array. When t > 0 and fi = v = 0 in equation (7), the 
state \F) can be written, in the spin representation, as 



\F) = 



iV2\nU---i) IfN = even, 
\v5|tm---t) HN = odd, 



where 



lt) = 



, 14) = 



(31) 



(32) 



are the two eigenstates of a z with eigenvalues 1 and — 1, respectively. 

It can be derived that yf = o\ and y B N = — i<J z N Pf, where the parity 
operator associated with the Z 2 symmetry of the system is given by 



Pf= Hal. 



(33) 



In the spin representation, the two degenerate ground states and 
I^r) can be written as 

I^HyflFHaflF) 

_ f V2\UU ■ ■ ■ l) IfN = even, (34) 

l^lTm-"T> IfN = odd, 

and 

\V R )=y B N \F) = -ia z N P f \F) 

-/ v / 2|iTiT---T> IfiV = even, (35) 
.*V2|mt"4) IfN = odd. 

Also, it is clear that 0F L | = 0. The two basis states |0) = d end \F) 
and 1 1) = dg nd |0) of the Majorana qubit are given, respectively, by 

\o) = \(vi+*&)\f) 

1 -(itm--a>+imt---T)) ifN= ev en, (36) 

= (IUU---T>-imT-"i» ifN=odd, 
|i> = ^(yf-^)|o> = ^(«?-^/)|o> 

:(IUU--a>-|mt"-t» IfN = even, (37) 

=(iTm-"T>+imT--i» ifN=odd, 



and 



These two basis states of the Majorana qubit are also two degenerate 
ground states of the system and have well-defined parities. 

When the parameters are identified at each qubit, we can rewrite 
equation (7) as 

N-l N 
n=l n=l 

Here v n = I p ® 0 Q -f)j , and f n = <D W /<D 0 H (n) / 2 > where 
=®( n ) J 0> 0 . The interqubit coupling reads 

or An,n+\) ±(n,n+l) 

PEj c cos [In ft -<Pc 



tn,n+l : 



i+2 i Scos(27 t / c ( "'" +1) -^"-" +1) )' 



(39) 
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where 



(«,«+!) . 



2Psm(2nf c {n > n + 1] 
l + 2^cos(27i/ c (n ' n+i: 



(40) 



and ^(n,n+ 1) = q,(n,n+ 1) j q q is the re( j UC ed fi^ applied to the coupler 

between qubits w and n + 1. 

The local perturbation (5/i„ cr* can only tend to drive the ground 
state (i.e., the Majorana-qubit state) |0) (| 1)) to an excited state of the 
system which has an energy level higher than the ground state. 
Similar to the case of the charge-qubit array, this is also owing to 
the protection of the Majorana-mode states and against 
the local perturbation 3/n n <f n . Indeed, the local perturbation 3\i n 
tends to drive |*F L ) (I^Pj*)) to an excited state with an energy differ- 
ence A from the ground state, where A = At for 1 < n < AT, and A = 
2t for n = 1 and N. Nevertheless, such a state transition is not 
permitted for a small perturbation 3/n n o x n . Therefore, in contrast 
to a single flux qubit, the local perturbation Sfi n o x n will not produce 
decoherence to the Majorana qubit. 

The local perturbation 3v n o z n randomly shifts the energy levels of 
the Majorana-mode states and to cause pure dephasing to 
these states. Also, it produces relaxation to the Majorana-qubit states 
|0) and 1 1). However, the coupler perturbation 3t n , n+ i o z n o z n+l yields 
pure dephasing to both the Majorana-mode states Q^i) and 1^)) 
and the Majorana-qubit states (|0) and |1)). Because the interqubit 
coupling is usually much smaller than I p Q>o, the coupler perturbation 
St n>n+ i a z n a z n+l should be much weaker than 3v n o z n . Therefore, in 
the case of a flux-qubit array, the dominant noise of the Majorana 
qubit is due to the perturbation 3v n a z n . In order to improve the 
quantum coherence of the Majorana qubit, one can suppress the 
fluctuations 3v n by reducing I p . This can be achieved by reducing 
the size of the Josephson junctions in each flux qubit because I p is 
proportional to the Josephson coupling energy Ej. Note that when 
reducing the size of the Josephson junctions to suppress flux noise, 
the charge noise can finally become important, due to the increasing 
charging energy. Thus, while the flux noise is suppressed, one can 
shunt a large capacitance to the Josephson junction, so as to suppress 
the charge noise as well. This method was proposed to increase the 
decoherence time of the flux qubit 42 and has been implemented in a 
recent experiment 43 . 

Manipulating and probing Majorana modes. The superconduct- 
ing-qubit arrays proposed above can be used to realize a tunable ID 
quantum Ising model on wire networks, similar to the semicon- 
ducting wire networks in Ref. 16, to demonstrate the non-Abelian 
statistics of Majorana fermions. In particular, braiding Majoranas 
can be implemented via a T-junction formed by two perpendicular 
wires 16 . Here we use four superconducting qubits, as the smallest size 
of the system, to form such a T-junction [see Figure 3(a)], where v = 
0 for all charge (flux) qubits. When the Jordan-Wigner tranfor- 
mation is performed, this T-junction of four qubits is described by 



H = t{ai—a[ s j {a\ + a 2 ) + 1 (ay - a\)j {a\ + a 2 ) 
+ 1 (a 2 — <4) ( a l + ~ M ^2a|ai — 1^ 

— [i(2a\,ai> — 1^ — ii(la\a 2 — 1^ — ^2a|a 3 — 1^ 

- i {wiY B i + vriyy + ^2 yl + wiv B ) , 

where qubits are numbered by starting from sites 1 and 1 ' and ending 
at site 3. 



For the Hamiltonian (41), when t = 0 for all pairs of adjoining 
qubits, the nearest-neighbor Majoranas related to the same site are 
coupled by fi, and the nearest-neighbor Majoranas related to two 
adjoining sites are decoupled. Then, Hamiltonian (41) is reduced to 



H = 



(42) 



In this case, the Majoranas are all paired in the whole T-junction 
region and no edge states occur. Starting from this phase, we adia- 
batically vary the parameters of superconducting qubits to have the 
horizontal array become an unpaired- Majarana region, i.e., adiabat- 
ically turn the parameter fi to zero for each qubit in the horizontal 
array and simultaneously switch on the interqubit coupling t for the 
horizontal array. Then, the Hamiltonian (42) becomes 



H = i(tf 1 y* + ty B 2 yl)-inriy B r . 



(43) 
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Figure 3 | Braiding two unpaired Majorana fermions. (a) T-junction 
formed by qubits 1, 2, 3 and 1', where each qubit is denoted by a 
rectangular box. Two Majoranas related to the same qubit (e.g., and yf ) 
can be paired by the parameter fi while two Majoranas related to adjoining 
qubits (e.g., and y£) can be paired by t. (b) Unpaired-Majorana region 
for the whole horizontal array, where an unpaired Majorana (denoted by a 
solid circle) is located at each end. (c) Adiabatically tuning fi to nonzero for 
qubit 1 and turning off f between qubits 1 and 2 drive the left-end Majorana 
mode (shown in orange) to the middle qubit. (d) Adiabatically tuning fi to 
zero for qubit 1 ' and turning on t between qubits 2 and 1 ' drive the original 
left-end Majorana to the bottom of the T-junction, (e) The right-end 
Majorana (in black) is driven to the middle qubit by adiabatically tuning /u 
to nonzero for qubit 3 and turning off t between qubits 2 and 3. (f) The 
original right-end Majorana is finally driven to the left end by adiabatically 
tuning /u to a sufficiently large value for qubit 2, tuning /u to zero for qubit 1 
and turning on t between qubits 1 and 2. (g) The Majorana at the bottom is 
driven to the middle qubit by adiabatically tuning /u to zero for qubit 2 and 
tuning fi to a sufficiently large value for qubit V . (h) Adiabatically turning 
off t between qubits 2 and 1 tuning /u to zero for qubit 3, and turning on t 
between qubits 2 and 3 finally drive the original left-end Majorana to the 
right end. This accomplishes the anti-clockwise braiding of two Majorana 
fermions. 
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This corresponds to the configuration of Majoranas in Figure 3(b), 
where a pair of isolated Majoranas emerge at the two ends of the 
horizontal array. Here adiabatic changes of the parameters with 
respect to the time z require that 

H\d^/dz\<^E 2 g , h\dt/dz\^E 2 g , (44) 

where E g is the energy gap between the first excited and ground states 
of the system at the time z. Generally, the system takes the super- 
position state of these two degenerate Majorana modes: 
|*F) = (ayf + fryf ) \F), where \F) is the state in which all eigenstates 
of the system with E < 0 are occupied. However, while reaching the 
state in Figure 3(b), if p for qubits 1, 2 and 3 are all adiabatically tuned 
to zero in the same manner and the interqubit coupling between 
qubits 1 and 2 is adiabatically switched on in the same way as that 
between qubits 2 and 3, then the left- and right-end Majoranas 
should occur with equal probabilities. Using this state 

|*F).= - (yf + e ie y%)\F) as the initial state, one can braid the left- 

and right-end Majoranas through the steps shown in Figures 3(c)- 
3(h) by adiabatically tuning the qubit parameters. For instance, by 
adiabatically switching on p for qubit 1 and turning off the coupling t 
between qubits 1 and 2, Hamiltonian (43) is changed to 

H=-i( W fyf + W f.yf.)+i^ ) (45) 

i.e., the configuration of Majoranas in Figure 3(b) is adiabatically 
converted to the configuration of Majoranas in Figure 3(c). 
Similarly, other steps shown in Figures 3(d)-3(h) can be achieved. 
This braiding of Majoranas following the steps from Figure 3(b) to 
3(h) corresponds to a unitary operator 416 which transforms yf to yf 
and yf to — yf . Therefore, the initial state I*?),- of the system is 

transferred to \ x ¥) f = - (yf -e i6 yf ) \F) after braiding the left- and 

right-end Majoranas. 

Finally, we focus on probing Majorana fermions. The initial state 

|*F) ; = - (yf — e^yf ) \F) given in Figure 3(b) is a ground state of the 

system with qubit 1 ' decoupled from the horizontal array of super- 
conducting qubits. When expressed in the basis states of qubits, this 
initial state can be written as I*?),- = | x J / 12 3>f ® where 

| X Fl23>^^l|0l0 2 0 3 )+^2|0i0 2 l 3 > 

+ ^ 3 |0 1 1 2 0 3 ) +A/ 4 |0 1 1 2 1 3 ) 

(46) 

+ ^5|li0 2 0 3 >+^ 6 |l 1 0 2 l 3 > 

+ ^7|lll20 3 )+^8|lil 2 l3>, 

and |*Fr), = in\0y) + £ /2 |li'}- Also, the final state 
\ x ¥)f= ~ (yf — e ld y\ )\F) is another degenerate ground state of the 
same system and can be expressed as |*F)y = ^123)/ ® l^i')^ where 
|**V)/ = <aI°i'> + ^l 1 !')' and |*F 123 )/has the same form as |*F 123 )„ 
but the Xn are replaced by Xp with / = 1 to 8. The states \ x ¥) i and \ x ¥)f 
can be distinguished using experimentally available state-tomo- 
graphy techniques for superconducting qubits (see, e.g., Refs. 
44,45), which involve reconstructing an unknown quantum state 
from a complete set of measurements of the system observables. 
For the initial state in Figure 3(b) and the final state in Figure 3(h), 
the qubit V is decoupled from the array consisting of the three 
coupled qubits 1, 2 and 3. Thus, only one-qubit tomography for qubit 
1' and three-qubit tomography for coupled qubits 1, 2 and 3 are 
required for distinguishing the initial and final states. 

Note that the initial and final states after braiding yf and yf can 

be written as |«F>,= \ (yf \F) +e i9 yf \F}) and |T) /= 1 (yf \F) - 
e t6 yf\F)), which have different relative phases between yf\F) and 



yf \F). Because the qubit V is decoupled from the horizontal array 
consisting of the three coupled qubits 1, 2 and 3, the states |*Fi')i and 
|*Fi')y are the same, in addition to a global phase between them. Thus, 
when performing quantum-state tomography, the different relative 
phases between yf\F) and y%\F) will give rise to the difference 
between 1%^}; and |¥i 23 )/. 

Experimentally, it is more complicated to use state-tomography 
techniques to determine the quantum state of three qubits other than 
two qubits. Therefore, we can consider the state | *F) . in Figure 3(c) as 
the initial state. This state is a ground state of the system with qubits 
1' and 1 decoupled from other qubits and can be decomposed as 
1*),.= 1^23),® l^r),®!^),-, where 

1^23} = An |0 2 0 3 ) + la |0 2 1 3 ) + Afl 1 1 2 0 3 ) + In 1 1 2 1 3 ) , 
l^r^&ilOrJ + kllr). l^i>« = 'Jfl|Oi)+iJ 0 |li>. ^ 

From the state in Figure 3(h), further proceeding with one step 
analogous to that from Figure 3(b) to Figure 3(c), we achieve 
the final state with the originally unpaired Majoranas yf and yf 
braided. This final state can also be decomposed as 

|*) / = |^23> / ®|^r) / ®|^i> / , where \¥ v ) f = ^\0 V ) + ^|l r ), 

l^i)/ = Vfi\0i) + ^|liX and |*F 23 )/has the same form as \ x ¥ 2 3)v 
but the Xi\ are replaced by Xp with / = 1 to 4. Similarly, the states |*F) . 
and |*F)^. can also be discriminated using state-tomography tech- 
niques. Here, because qubits 1 and 1 ' are decoupled from the two 
coupled qubits 2 and 3, only two-qubit tomography for coupled 
qubits 2 and 3 as well as one-qubit tomography for qubits 1 and 1 ' 
are needed for distinguishing the initial and final states. 

Experimentally, in addition to one-qubit tomography, two-qubit 
tomography is also implementable for superconducting qubits (see, 
e.g., Refs. 44,45). Thus, it is feasible to measure and |^}^ 

because the two-qubit tomography can be used to determine 1^23)1 
and 1^3)/. Moreover, quantum-state tomography has been per- 
formed on three 46 or even five 47 superconducting qubits, so it also 
becomes feasible to measure and l^/by determining | ¥123)1 
and | x F 123 )y via quantum-state tomography. This is important here 
since information might be lost by only performing two-qubit tomo- 
graphy, particularly in the case of poor gate fidelity or decoherence. 

Discussion 

When fabricating superconducting circuits, parameter variations 
unavoidably occur, as in any solid-state system. For the charge-qubit 

array, v = 0 can be achieved by having <& q = ~®o> irrespective of the 

parameter variations. Also, /n can be tuned, via the gate voltage V g , to 
the required value, even if E ch varies for different qubits. For varying 
E J0 among qubits, the interqubit couplings also vary [see equation 
(2)]. One can replace the large Josephson junction by a dc SQUID 
and tune the SQUID, i.e., the effective E Jc , to obtain the desired 
value t for the interqubit coupling. For a symmetric SQUID with 

Josephson coupling energy the effective E Jc is given by 

£; c (O c ) =2£r^cos(7cO c /(Do)» where <J> C is the magnetic flux in the 
SQUID loop. In equation (2), E Jc is now replaced by E Jc (<b c ); in both 

equation (2) and v, <S> q is replaced by = <S> q + - <$> c . Therefore, the 

tunability of t(Q> c /5>'qj with v = 0 can be implemented by changing 
both (J> c and Q> q . 

As for the flux-qubit array, v = 0 can be achieved by having/ = - . 

Also, pi can be tuned to the given value by changing the flux/ s applied 
to the SQUID in each qubit. Moreover, even if the parameters of 
couplers vary, one can tune the flux/ c in each coupler to achieve the 
required value of t for the interqubit coupling [see equation (8)]. 
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Furthermore, note that even if v = 0 cannot be experimentally 
reached very accurately, our proposal still works for the states in the 
unpaired-Majorana region of the system if |v/t|«Cl. Experimentally, 
this can be achieved by designing a relatively strong interqubit coup- 
ling for the qubit array. For instance, because I p ~ 2nEj/<$) 0 , one has 
1 



\v/t\~[2n(l + 2P)Ej/PEjc 



When jff = 0.1, E /c = 5E P and/ * 



[0.499, 0.501], |v/f|~0.015<l. 

In conclusion, we propose superconducting circuits to construct 
two superconducting- qubit arrays where Majorana modes can occur. 
The unpaired zero-energy Majorana modes, which emerge at the left 
and right ends of the chain in the Majorana-fermion representation, 
can be used to encode a qubit called the Majorana qubit. Also, we 
express this Majorana qubit in the spin representation and show its 
advantage, over a single superconducting qubit, for quantum coher- 
ence. Moreover, we suggest using four superconducting qubits as the 
smallest circuit to demonstrate the braiding of Majorana modes, and 
show how to distinguish the states before and after braiding 
Majorana modes. These superconducting- qubit arrays can, in prin- 
ciple, be extended to wire networks, similar to the semiconducting 
wire networks in Ref. 16, to demonstrate the non-Abelian statistics of 
Majorana modes. 
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